Constraints on a scalar-pseudoscalar Higgs mixing at future e + e colliders: an update 
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We perform an update of our previous analysis of the constraints on possible deviations of Hbb 
coupling parametrized as —(0 + 1756), arising from a scalar-pseudoscalar mixing, where the process 
e + e~ — > bbvv was used. In this paper we include a complete simulation of the process e + e~ — ► 
bbe + e~ and combine these results to obtain tighter bounds on the deviations of the parameters a 
and b from their Standard Model values that could be measured at the Next Linear Collider. 
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I. INTRODUCTION 



The origin of fermion masses and mixings is an out- 
standing open problem in particle physics. In the Stan- 
dard Model (SM), the Higgs mechanism is responsible for 
the electroweak symmetry breaking and mass generation 
via ad-hoc Yukawa couplings. There are reasons to be- 
lieve that the SM is not the final model and a complete 
study of the coupling of the lightest boson, which we will 
call the Higgs boson, to fermions can provide hints on 
new physics beyond the SM. 

In a recent letter [1], we performed a realistic simu- 
lation of the process e + e~ — > bbvv, where v can be an 
electron, muon or tau neutrino, in the environment of a 
future Linear Collider with a center-of-mass energy of 
y/s = 500 GeV with an accumulated luminosity of 1 
ab _1 , based on the TESLA design [2]. In particular, 
we noticed that weak gauge boson fusion is the domi- 
nant contribution to the subset of diagrams containing 
the Higgs boson for M H < 180 GeV at y/s > 500 GeV 
and hence this process is sensitive to the Higgs boson cou- 
plings to &-quarks. We have used a similar technique to 
investigated the possibility of detecting deviations from 
the SM in the Higgs couplings to r-leptons at future e + e~ 
colliders [3], which can be improved by using r spin cor- 
relations [4]. 

In this Brief Report, we update the result [1] by com- 
bining it with the results of a detailed analysis of the pro- 
cess e + e~ — > bbe + e~ . This process has a very clean final 
state, with no missing energy and easy reconstruction, 
which can compensate for its smaller rate. It certainly 
must be included in a global analysis of the Hbb vertex. 

In order to perform our analysis, we will assume that 
the Higgs boson has already been discovered at the Large 
Hadron Collider and concentrate on the determination of 
its coupling to &-quarks. In extensions of the SM with ex- 
tra scalars and pseudoscalars, the lightest spin-0 particle 
can be an admixture of states without a definite parity. 



Hence, we parametrize the general Hbb coupling as: 



(a + 175 b) 



(1) 



where v = 246 GeV, m b is the 6-quark mass and a = 1, 
b = in the SM. In CP-violating extensions of the SM, 
deviations of these parameters may be generated at tree- 
level [5]. 

We will present results considering a and b as indepen- 
dent parameters and also for the cases of fixed a = 1 , free 
b and fixed b = 0, free a. There is a region of insensivity 
around circles in the a — b plane since we can't at this 
level of analysis disentangle the effects of a and b. 

In the SM the e + e~ — > bbe + e~ process is determined 
by 50 Feynman diagrams. Only 2 of these (Higgs ra- 
diative production, Z* — > > bb)Z, and vector boson 
fusion, Z*Z* — > H(— > bb)) can be considered as signal, 
where deviations from SM Higgs couplings to 6-quarks 
show up. The remaining diagrams are not changed by 
new physics in the Higgs sector. 

The total SM cross section for the process e + e~ — > 
bbe + e~ is approximately 43 fb for Mh = 120 GeV at 
yfs = 500 GeV, with cuts in the scattering angle be- 
tween initial and final electrons, | cos 6 ee \ < 0.9962, and 
in the final electron-positron invariant mass, M e + e - > 2 
GeV. The angular cut avoids the region of ±5° around 
the pipeline, used by the beam pipe itself and by the 
Small Angle Monitor (SAM) for luminosity measure- 
ments. The invariant mass cut insures that the elec- 
trons and positrons will be detected in the lead glass 
calorimeter [2]. Backgrounds from e + e~ — > e + e~ ZZ — > 
e + e~bbi>i>, e + e~ — > ZZZ — > e + e~bbvv, etc, can be re- 
duced to the 0.1 fb level by appropriate cuts [6]. 

In our analysis we will assume SM couplings of the 
Higgs to the electroweak gauge bosons. A full analysis 
taking into account non-SM HWW and HZZ couplings 
would introduce many more parameters and is beyond 
the scope of this paper. 
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II. ANALYSIS AND RESULTS 



We performed our Monte Carlo simulation by gener- 
ating observables represented as series in the a and b 
couplings multiplied by kinematical factors: 



= A + a ■ A x + a 2 • A 2 + ab ■ A 3 + 
b ■ A 4 + b 2 ■ A 5 . . . 



(2) 



FIG. 2. Contribution of Higgs diagrams to the differential 
cos# e (, distribution in the e + e~ — * bbe + e~~ process for Stan- 
dard Model (o = 1,6 = 0, black dots) and a — 1.0,6 = 0.5 
(crossed dots). 



where O is any observable and the Ai terms are purely 
kinematical structures which do not contain any a and b 
dependence and results from the amplitude squaring and 
phase space integration. In this way we only simulate 
the process once for each observable. In our case, terms 
linear in b parameter vanish and hence A3 = A4 = 0. 
We have studies the following observables: transverse b- 
quark momentum pTb, bb invariant mass M bb , cos9 eb , 
where 9 eb is the scattering angle between the &-jet and 
initial beam directions, and the T correlation, defined by 

T = JJSrPei ■ (PbXPb)- 

The event sample reproducing the expected statistics 
at TESLA was generated using our Monte Carlo pack- 
age while the detector response was simulated with the 
code SIMDET version 3.01 [7]. We assume an efficiency for 
6— jet pair reconstruction of e bb — 56 %, which is based 
on the 6-tag algorithm, as assumed in ref. [6]. In our 
simulations we used Mh = 120 GeV. 

In Figure 1 we show, for comparison purposes, the dif- 
ferential distributions in the invariant bb mass and cos 9 ebl 
for the total SM contribution and for the Higgs contri- 
bution only (including interference with SM). Notice the 
Z-boson peak and the smaller Higgs peak in the M bb 
distribution, and how the Higgs contribution dominates 
around its peak. 
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FIG. 1. The differential cross sections for cos9 eb and M bb 
for the process e + e~ — > 66e + e~. Solid circles are the full SM 
result while crossed circles are the contribution from Higgs 
only (including interference effects). 



In order to demonstrate the effect of different values of 
the parameters a and b, we show in Figure 2 the cos 9 e b 
distribution arising from the Higgs contribution for the 
SM (a = 1, b = 0) compared to the case with a = 1, 
b = 0.5. We see that the shapes are very similar, as 
expected, but the levels can be noticeable different. 

Another important aspect is the assumption about the 
detector performance and possible sources of the system- 
atic uncertaintes. We include the anticipated systematic 
errors of 0.5% in the luminosity measurement, 1% in the 
acceptance determination, 1% in the branching ratios, 
and 1% in the background substraction, and assume the 
Gaussian nature of the sytematics [2] . To place bounds 
on the Hbb couplings, we use a standard x 2 -criterion 
to analyse the events. After various kinematical distri- 
butions were examined, we found that the most strict 
bounds are achieved from cos 9 eb distribution by dividing 
the distribution event samples into 10 bins. The experi- 
mental error Aa^ xp for the i th bin is given by: 



Act, 



exp 



SM 



\J <>syst + <>st. 



where 



^<jf M e bb jCdt 



(3) 



(4) 



and Sg yst is the sum in quadrature of the systematic un- 
certainties mentioned above. 

Below we present our final results for a TESLA-like 
environment [2] with a center-of-mass energy of 500 GeV 
and for Mh = 120 GeV. We investigated three possible 
scenarios for the luminosities: 100 fb _1 , 1 ab -1 and 10 
ab" 1 . 

The bounds that can be obtained at 95% confidence 
level from e+e~ — > bbe + e~ process on the Aa = a — 1 
and b parameters are: 

-0.09 < Aa < 0.08 for C = 100 fb" 1 ; (5) 

-0.056 < Aa < 0.055 for £ = 1 ab" 1 ; 

-0.05 < Aa < 0.05 for C = 10 ab~\ 

for the case of b = and free Aa and 

-0.42 < b < 0.42 for C = 100 fb" 1 ; (6) 

-0.32 < b < 0.32 for C = 1 ab" 1 ; 

-0.3 < b < 0.3 for C = 10 ab~\ 
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for the case of Aa — and free b. 

We will combine these limits with our previous bounds 
from the e + e~ — > bbvv process according to the following 
procedure. Given two bounds on the same quantitity X, 
say \X\ < a and \X\ < c 2 , the combined bound will be 



\X\ < c 3 where c 3 = ^ (ci^i) 2 + (c 2 o> 2 ) 2 , 
and the statistical weights u)\ and w 2 are given by 
N\ , A 2 



JVi + iV 2 



and c^2 = 



(7) 



(8) 



where N\ and A^ 2 are the number of events for each pro- 
cess. 

With these prescription we obtain the improved 
bounds: 



-0.037 < Aa < 0.035 for C = 100 fb" 1 ; 
-0.024 < Aa < 0.024 for C = 1 ab" 1 ; 
-0.022 < Aa < 0.022 for C = 10 ab"\ 



(9) 



are comparable to the study performed in [2], where a 
global fit analysis for C = 500 fbr 1 and ^ = 500 GeV 
has resulted in a relative accuracy of 2.2% in the gnbb 
Yukawa coupling. For comparison, the top quark Yukawa 
coupling can be determined with a statistical accuracy of 
16% at the LHC for M H = 130 GeV [8]. 

In our analysis we can not disentangle the contribu- 
tions from deviations in the a and b parameters. How- 
ever, notice that the partial width T H ^ bb - is proportional 
to (a 2 + b 2 ). Our distributions have a different depen- 
dence 



da 
dO 



= A Q + a ■ Ai + a 2 • A 2 + b 2 ■ A 3 



•bb 



Therefore, if an independent measurement of 
obtained (for instance, from on-mass-shcll Higgs pro- 
duction in Higgs-strahlung or in a muonic collider), one 
would be able to separate the a and 6 contributions and 
obtain an explicit indication of CP violation in the Higgs 
sector. 



for the case of b = and free Aa and 

-0.24 <b< 0.24 for C = 100 fbr 1 ; 
-0.20 < b < 0.20 for C = 1 ab" 1 ; 
-0.19 < b < 0.19 for C = 10 ab" 1 , 



(10) 



for the case of Aa = and free 6. 

These results can be roughly scaled for moderate vari- 
ations in the Higgs boson mass around 120 GeV by mul- 
tiplying the bounds by a factor (M H /120 GeV) 2 . 
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III. CONCLUSIONS 

We have performed an update in our previous con- 
straints on deviations of the Hbb coupling from its SM 
value by including a complete analysis of the sensitivity 
due to the process e+e~ — ► bbe + e~ at the next genera- 
tion of linear colliders. These deviations are predicted by 
many extensions of the Standard Model. We improved 
our previous bounds by roughly 10%. 

We showed that future e + e~ linear collider experi- 
ments will be able to probe deviations of Hbb coupling. 
The weak gauge boson fusion process is instrumental for 
achieving such a precision. For a TESLA-like environ- 
ment, we are able to constrain the couplings at the level 
of a few percent for the a parameter (for fixed b) and tens 
of percent for the b parameter (for fixed a). These results 
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